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Summary We present a theory of homogeneous volatility bridge estimators for log- 
f-H , price stochastic processes. The main tool of our theory is the parsimonious encoding 

. of the information contained in the open, high and low prices of incomplete bridge, 

• ' corresponding to given log-price stochastic process, and in its close value, for a given 

C ' time interval. The efficiency of the new proposed estimators is favorably compared with 

, that of the Garman-Klass and Parkinson estimators. 

Keywords: volatility, variance, estimators, efficiency, Wiener processes, homoge- 
neous functions 

> ■ 1. INTRODUCTION 

Volatility, defined as the standard deviation of the increments of the log-price over a 
\^ . specific time interval, is a universally used risk indicator. With the growing availabil- 

ity of high-frequency tick-by-tick price time series, a number of new efficient volatility 
. estimators have been developed (see, for instance, Yang and Zhang (2000), Corsi et al. 

(2001), Andersen et al. (2003), A'l't-Sahafia (2005), Zhang et al. (2005)). However, for 
most applications involving risk assessment and management of investment portfolios, it 
is the common practice to use time series of prices recorded at fixed time intervals, such 
as 1 minute, 5 minutes, 30 minutes, 1 hour, 1 day, 1 week, 1 month and so on. For such 
k> \ time series, four prices are actually recorded, called the open-high-low-close (OHLC) of 

. the price for each time interval. 

' Our purpose is to provide new tools to exploit systematically the OHLC to improve 

volatility estimators, compared with techniques using only the close price time series. 
It is intuitively appealing that close-minus-open, high-minus-open and low-minus-open 
should provide significant information on the variability of the price, that should help 
improve the volatility estimators. We present here a comprehensive theory of homoge- 
neous volatility bridge estimators for arbitrary stochastic processes, that fully exploit the 
OHLC prices. For this, we have started to develop the theory of the most efficient point- 
wise homogeneous OHLC volatility estimators, valid for any price processes (Saichev et 
al., 2009). The main tool of our theory is the parsimonious encoding of all the informa- 
tion contained in the mentioned OHLC prices for a given time interval in the form of 
general "diagrams" associated with the joint distributions of the high-minus-open, low- 
minus-open and close-minus-open values. The diagrams can be tailored to yield the most 
efficient estimators associated to any statistical properties of the underlying log-price 
stochastic process. 

The present work extends and generalizes (Saichev et al., 2009) by developing most 
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efficient OHLC bridge estimators. We find that the new OHLC bridge estimators are 
significantly more efficient than the OHLC estimators obtained from the untransformed 
(or "unbridged" ) process. For Wiener and similar processes, this can be intuitively un- 
derstand as follows. It is well-known that the high and low values of a Wiener process 
are most probably found in the neighborhood of the edges of the observation interval. 
In contrast, by construction of the bridge, its high and low values are in general distant 
from the edges. As a result, the high and low of a bridge incorporate significantly more 
information about the variability of the original stochastic process than its own high and 
low values. This is the motivation for us to extend the theory of (Saichev et al., 2009) 
for bridges and to provide explicit analytical expressions for the most efficient point-wise 
volatility bridge estimators, based on the analytical expression of the joint distribution 
of its high-minus-open, low-minus-opcn and closc-minus-open values. 

Our work also improves on the following papers as follows. Garman and Klass (G&K) 
(1980) introduced a quadratic estimator for the variance of the Wiener process for the 
log-price, which has rather low variance. Parkinson (PARK) (1980) proposed a simple 
quadratic variance estimator proportional to {H — i)^, which is using only a part of 
the information available from OHLC prices. Rogers and SatchcU (R&S) (1991,1994) in- 
troduced another quadratic estimator for the variance of the Wiener process with drift, 
which is unbiased for all drifts. Both G&K and R&S estimators are focused on the 
variance, and do not present estimators for the volatility, which is of obvious interest 
for financial applications. Yang and Zhang (2000) produced an unbiased and efficient 
quadratic variance estimator, taking into account the OHLC of log-prices for n > 1 con- 
secutive days. Their main novelty is to take into account the possible existence of jumps 
(or gaps) of prices from yesterday's close till today's open prices. Their minimization of 
the variance of their estimators requires the estimation of expectations of a quadratic 
form of the OHLC which they only partly achieve due to the lack of knowledge of the 
full joint distribution, which we offer in this paper. Chan and Lien (2003) compared the 
empirical effectiveness of four estimators, the PARK, the G&K and R&S ones, and the 
naive excursion range H — L estimator. From the perspective offered by these previous 
works, the present paper can be viewed as providing their full underpinning theory, since 
we are able to express efficient estimators in the presence of arbitrary constraints from 
the explicit knowledge of the joint distribution of the OHLC log-prices. 

The paper is organized as follows. Section 2 describes the properties of the stochastic 
processes for which our theory of most efficient homogeneous volatility bridge estimators 
is developed. Section 3 derives the general expressions for the most efficient homogeneous 
volatility OHLC bridge estimators. Section 4 provides a detailed analytical description of 
the statistical properties of incomplete bridges of Wiener process with drift, describing 
log-price dynamics. Section 5 compares the efficiency of our derived most efficient ho- 
mogeneous bridge estimators and the efficiency of the generalized G&K bridge estimator 
and of the normalized PARK estimator. Section 6 tests our results using synthetic time 
series generated by numerical simulations, which mimic the tick-by-tick nature of real 
log-price processes. Section 7 concludes. 

2. HOMOGENEOUS VOLATILITY BRIDGE ESTIMATORS 

The main goal of this paper is to construct efficient bridge estimators using the open, 
high, low, close (OHLC) prices for the variance and the volatility of some asset log-price 
process A{t). 
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2.1. Volatility of order A 

The conventional definition of the volatility V{tQ, Tq) of a stochastic process A{t) at time 
to and time scale Tq is the standard deviation of its increment 

A(to,To) = ^(to + To)-A(to) 

within the time interval t G (<o, + ^o): 

y(to,To) = v/Var[A(to,To)]. 

The time scale Tq can be for instance 5 minutes, 1 day or 1 year, corresponding respec- 
tively to intraday, daily or yearly volatility. 

Since different measures of the variability of log-price processes are used in the litera- 
ture, it is convenient to define a generalized volatility of order A as follows. 



Definition 2.1 The volatility of order A of the stochastic process A{t) is the power A of 
the conventional volatility 

Fa (to, To) := V\tQ,TQ) = (Var [A(to, To)])^/' . 



Remark 2.1 For A = 1, the volatility of order A coincides with the conventional volatil- 
ity, while, for A = 2, V2(to,2o) is the variance of the increment A(io,7o)- Most known 
estimators, for instance the R&S, G&K and PARK ones, are variance estimators. Intro- 
ducing the volatility of order A gives us the possibility later on to compare the differences 
and relations between the volatility and variance estimators. 



2.2. Wiener process model of log-price increments 



We will analyze the properties of the estimators of the volatility of order A for the Wiener 
process with drift, posing without loss of generality to = and A{Q) = 0. This implies 
that 

A{t) -.^ tit + <jW{t), (2.1) 

where fi is the drift of the log-price process A{t) and a is its standard deviation at t = 1, 
while W{t) is the standard Wiener process with zero drift and variance E[VF^(t)] = t. 
The self-similar properties of the Wiener process allow us to choose the time scale by 
To = 1 without loss of generality, so that the volatility of order A is simply equal to the 
standard deviation a raised to the power A: 

Vx:^Vx{tQ^O,TQ = l)=a\ 

We analyze below the volatility estimators based on the high, low and close values of the 
incomplete bridge of the Wiener process with drift A{t) defined by (2.1). 



Definition 2.2 The stochastic process 



B{t, K, T) A{t) - K—A{T) ^ ^(1 - K)t + a 



W{t) - n-W{T) 



(2.2) 



where k is arbitrary constant, is called the incomplete bridge of the original stochastic 
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process A{t). For k = 1, the incomplete bridge is nothing but the standard (complete) 
bridge 



B{t,T) := a 



Using the self-similar properties of the Wiener process, one can rewrite (2.1), (2.2) in 
the form 



(2.3) 



A{t) - aVrX (1,7), B{t, T) = aVrr 7^ , 

where y(t, k, 7) is the incomplete bridge 

y(i,K,7) :==X(i,7)-KtX(l,7) (2.4) 
of the Wiener process with drift 

Xit,j):^jt + W{t), te(0,l), (2.5) 
and the auxiliary parameter 

7 = -Vt (2.6) 
a 

plays the role of a I''' standardized moment (or inverse coefficient of variation) of the 
distribution of increments of the process A{t) over the time interval T. Figure 1 shows a 
realization of the Wiener process W{t) and its complete bridge. The high and low values 
of the Wiener process and of its bridge are in general drastically different. 

Remark 2.2 In financial markets applications, both the drift ^ and the standard devia- 
tion a are unknown. Thus, the value of the parameter 7 is unknown as well. Nevertheless, 
for the convenience of our analysis, we will suppose in the following derivations that the 
value of parameter 7 is given. One can take into account the indeterminateness of the 
parameter 7 by exploring in detail the dependence as a function of 7 of the bias and of 
the efficiency of the OHLC volatility bridge estimators, following the analysis performed 
by Saichev et al. (2009) for k = 0. 



2.3. Homogeneous volatility bridge estimators 

Definition 2.3 A volatility estimator V\ is called an homogeneous OHLC volatility 
bridge estimator of order A if it has the form 

^x = j^h^{H,L,C), (2.7) 

where h\ is a homogeneous function of order A, the random values H and L, are the 
high and low of the incomplete bridge B{t, k, T) defined by (2.2) within the observation 
interval (0, T), 

i? := sup B{t,K,T), Z ~ inf B{t,K,T), 
te(o,T) te(o,T) 

and 

C := A{T) ^ fiT + aW{T) 
is the close value of the original stochastic process for the log-price A{t). 
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Figure 1. A realization of the Wiener process W{t) and complete bridge W{t) — tW{l) 



A remarkable property of homogeneous estimators defined by (2.7) is that, for a given 
7, their statistical properties do not depend on the duration T of the observation interval. 
Mathematically, this fact is expressed by the following theorem. 



Theorem 2.1 The estimator defined by (2.7) is equal to 

ax (T^hx{H,L,C), 



(2.8) 



where H and L are the high and low values of the incomplete bridge Y{t, n, 7) defined by 
expression (2.4) 

H := sup r(i,K,7), L :^ inf Y{t,K,-f), (2.9) 
te(o,i) *e(o,i) 

while C :— X{l,j) is the close value of the Wiener process X{t,'y) with drift, defined by 
(2.5). 

Proof. Substituting the right-hand-side of the equalities of (2.3) into the right-hand- 
side of expression (2.7) and using the homogeneity of the function h\, we obtain equality 
(2.8). ■ 



Definition 2.4 We refer to the function 

ex^hx{H,L,C) 



(2.10) 



as the canonical OHLC volatility bridge estimator of order A. Using this definition, one 
can rewrite expression (2.8) in the form 



(TA = o- ex- 



(2.11) 



Remark 2.3 The statistical properties of the canonical estimators (2.10) depend on the 
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standard deviation a only through the parameter 7 defined in (2.6), which we assume in 
the following derivations to be known. 

3. MOST EFFICIENT HOMOGENEOUS BRIDGE ESTIMATORS 
3.1. Diagrams of homogeneous bridge estimators 

It results immediately from expressions (2.11), (2.10) that the homogeneous estimator 
given by (2.7) is unbiased, if the expected value of the corresponding canonical estimator 
given by (2.10) is equal to unity: 

nex]^nhx{H,L,C)]^l. (3.1) 



Definition 3.1 The homogeneous volatility bridge estimator of order A given by (2.7) 
is called the most efficient one, for a given value 70 of the parameter 7 and for a fixed 
value of the parameter k, if, for 7 = 70 and fixed values of n and A, the equality (3.1) 
holds while the variance of the corresponding canonical estimator achieves the minimal 
value among the variances of all canonical estimators of given order A and for the same 
parameters 7 = 70 and k. 

In this section, we provide the explicit expressions of the most efficient homogeneous 
volatility estimators. For this, it is convenient to use a change of variables from the ran- 
dom variables {-ff, L, C} to their corresponding spherical (geographic) coordinate vari- 
ables {i?,e,$}: 

iJ = i?cosecos$, L = i?cosesin$, C = i?sine. (3.2) 
Inversely, we have 

R= JH^+L^ + C^, (=) = arrtan | , ^ ) , $ = arctan( — ). (3.3) 

Substituting (3.2) into (2.10) and taking into account the homogeneity of the function 
h\, we obtain 

eA = i?^ VA(e,$), (3.4) 

where 

ipx{0, (f) = /lA (cos 6* cos (/), COS 9 sine/;, sin0). (3.5) 



Definition 3.2 The function ip\{9,(p) defined by expression (3.5) is called the diagram 
of the canonical estimator of order A. 



Remark 3.1 The spherical coordinate system is intrinsic to homogeneous estimators, 
allowing us to split them into a known power function R and an arbitrary function 
of the variables O and $ (see Eq. (3.4)). The spherical coordinate system reduces the 
search of efficient OHLC estimators from three-dimensional functions to the appropriate 
two-dimensional function ip\{9, 0). 
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3.2. Domain of possible {0,$} values 

Below, wc will need the domain of existence for the values of the random variables 
{R, 9, $} defined by (3.3). First, it is obvious that R G (0, oo). The domain S of the pos- 
sible values of the two other random variables {O, $} depends on the interplay between 
the random high H and low L of the incomplete bridge given by (2.4), and the close 
value C of the Wiener process with drift defined by (2.5). It will be clear below that S 
depends on the parameter n. Thus, we denote it by S^. We use the same notation 5^ for 
the domain of the arguments {0, 0} of the diagram ip\{0, (f) defined by (3.5): {6*, (j)} € S^. 
Since iJ ^ and L ^ 0, in view of (3.2), we have 

L TT 

tan$ = — e (-00,0] ^ --s^$<0. 
H 2 

In turn, as it seen from (2.4) and (2.5), the values {iJ, satisfy to the inequalities 

L ^ (1 - k)C ^ H or, using (3.2), 

• * /sini'N^ / cos$ 

sin<P ^ (1 — Kj tanB ^ cos$ ^ arctan ^ B ^ arctan 



1 / sin \ ^ / cos $ \ TT I , , 

5k = <^ arctan 6 arctan , -- s$ $ < ^ . (3.6) 



Thus 

..._.Lr.taJ , 

1- kJ \1- kJ 2 

3.3. Most efficient OHLC homogeneous bridge estimators 

Let us denote the joint probability density function (pdf ) of the random variables {H, L, C} 
by Q{h, c; K,j). Then, the expected value of the canonical estimator defined by (3.4) 
is equal to 

E[ex\K,-f]=MxiK,^):^ JJijx{&^<l^)9x{0,^;K,j)cos9d9d^, (3.7) 

where 

/•oo 

g\{6,(l);K,"f) ~ / p^'^^Q{pcos6 coscj), pcos6 sin(j3, psinO; K,^)dp. (3.8) 
Jo 

Accordingly, at 7 = 70 and given k, one can represent the diagram of any unbiased, 
homogeneous estimator in the form 

ipx{0,(l);K,jo) ^ n(a a\ ^T^l^ — ^ '^IKJl' (^■^) 

jj G{d,(p) gA(c', 70) cos 6'd6'd0 

where G(0, 0) is an arbitrary function. The following theorem gives the expression for 
the diagram (3.9) corresponding to the most efficient (for any given k and 7 = 70) 
homogeneous estimator of order A. 

Theorem 3.1 The diagram of the most efficient (for a given k andj — "fo) homogeneous 
bridge estimator of order X is equal to 

^™e,A(0,0;^,7o) = ^4^7^^^' (3.10) 
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where 



Gx{e,(f>;K,j) = — -, £;,K,7 = // -— -cos0d0d<p. 3.11) 

g2\{9,(l);K,j) J J 52a(6',0;k,7) 

The proof is given in Appendix A.l. 

4. STATISTICAL DESCRIPTION OF INCOMPLETE BRIDGES 
4-1. Identical in law Wiener process 

In order to get the most efBcicnt homogeneous OHLC bridge estimator, we need the pdf 
Q{h,£,c; Kjj) of the high and low of the incomplete bridge Y{t,K,'y) defined by (2.4) 
and the close value C of the underlying process X(t,j) defined by (2.5). Before giving 
the explicit solution, it is useful to discuss their general statistical properties. 

Theorem 4.1 The incomplete bridge Y(t,K,j) given by (2.4) is identical in law to the 
diffusion process 

y{t,K,-f):=-fil~K)t + Wit,K), (4.1) 

where 

Wit, .) := (1 - t + (1 - .rt)W . (4.2) 

Proof. After substitution (2.5) into (2.4). we obtain 

Yit,K,-f)^^il-K.)t + Q{t,K), (4.3) 

where 

n{t,K) W{t) - K.tW{l). 

One can easily verify that Q{t,K) is a Gaussian process with zero mean and covariance 
given by 

i:[n{t,,K)n{t2,K)] = {t,At2)-[i-{i-fif]tit2, o^ti,t2^i. (4.4) 

Direct calculations show that the Gaussian process W(<, k) defined by (4.2) is also char- 
acterized by a zero mean and the same covariance (4.4). This implies that the incomplete 
bridge Y{t, k, 7) given by (4.3) is identical in law to the diffusion process y{t, k, 7) defined 
in (4.1). ■ 

4-2. Change of time 

Henceforth, for the analysis of the statistical properties of the incomplete bridge Y{t, k, 7) 
defined by (2.4), we will use the equivalence in law stated in theorem 4.1, which allows 
us to to replace the incomplete bridge by the diffusion process y{t, k, 7) defined by (4.1). 
As will be clear below, it is convenient to explore the extremal properties of the diffusion 
process y(t, k, 7) using the change of time 

^ = <'^'') ■■= i-t+n-nrt- 
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Inversely, 

t — t(T, k) := 



T+ (1 - k)2(1 - r)' 

The function T{t, k) maps the mterval t £ (0, 1) onto the same interval r G (0, 1). 
Let us introduce the auxiliary stochastic process 

Z{T,K,j):^yit{T,K),K,-f). (4.5) 

Using relations (4.1), (4.2), and self-similar properties of Wiener process, rewrite Z(r, k, 7) 
in the form 

^(^:^,7)= ^br + W{T)]. (4.6) 

T + (1 - k)"^(1 - r) 

Below we assume, for definiteness, k < 1. 

It follows from the construction (4.5) of the stochastic process Z{t, k, 7) and from the 

equality (4.6) that the following inequalities are equivalent 

L s;,y{t,K,j) H, ^ a + ar s^W{t) s^b + I3t, i, re (0,1), (4.7) 

where 

a = (1 - k)L, 6 = (1 - k)H, a = ^ —L 13= ^ —H - 7. (4.8) 

Additionally, the close value C = ^(1,7) of the stochastic process X{t,j) (2.5) is 
tied to the close value of the incomplete bridge Y{t, k,j) given by (2.4) by the equality 
y(l,K,7) = (1 — k)C. In turn, it follows from the identity in law of the stochastic 
processes Y{t,K,'f) and y{t,K,'-f) and from relations (4.5), (4.6) that one may replace 
y(l,Ac,7)by 

Z(l,K,7) = (l-«)[7 + I^(l)]. 

Thus, one obtains 

W/(l) = C-7. (4.9) 
4-3. Diffusion equation 

Let us define the probability 

/(/i, e, c; K, 7)dc := Pr{C e{c,c + dc)ne^ Y{t, k, 7) /i; t e (0, 1)}. 

Then, the joint pdf of the high and low values {H, L} (2.9) of the incomplete bridge 
Y{t, K, 7), and of the close value C of the original process X{t, 7) (2.5), is equal to 

Q{h,£,c;K.,j) = -^j-^^ , (4.10) 

£ h 

h> h^, i<£+, s$ c ^ , = V (1 - k)c, e+=OA{l- k)c. 

1 ~ K 1 — K 

From the relations of the previous subsection 4.2, one can express the function /(/i, £, c; k, 7) 
via the auxiliary function ip{u!;T) 

T)duj := Pr{VK(r) e {uj,uj + dui) n a + ar' s$ W{t') 5 + /3r'; r' G (0, r)}, 

according to 

f{h, £, c; K, 7) = if{c - 7; 1, a, 6, a, /3). (4.11) 
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The theory of Wiener processes iniphes that the auxihary function (^(o;; r) is the solu- 
tion of the diffusion equation 

(4 12) 

with initial condition 

Lp{oj-T ^Q) ^ 5{uj) (4.13) 
and absorbing boundary conditions 

(p{uj ^ a + aT]T) ^Q, f{uj ^ b + (3t;t) = 0, r > 0, (4.14) 

which account for the inequalities (4.7). 

Below, we solve this initial-boundary problem (4.12), (4.13), (4.14) and determine the 
joint pdf of the high and low values of the incomplete bridge Y{t, n,^) and of the close 
value of the Wiener process X{t,^) with drift, using the following relation that derives 
from (4.10) and (4.11): 

Q{h,e,c;K,-f) = ^j-^^ . (4.15) 

4-.4- Useful properties of the solutions of diffusion equations 



Before solving explicitly the initial-boundary problem (4.12), (4.13), (4.14), it is useful 
to present some general properties of its solutions. Firstly, if (/3(w; r) is some solution 
of diffusion equation (4.12), then Aipiuj + a;T), where a and A are arbitrary constants, 
is also a solution. Such relation tying together different solutions of the same diffusion 
equation can be written as 

(/3(w;t) < — > ALp{uo + a]T). (4.16) 

In order to solve the initial-boundary problem (4.12), (4.13), (4.14), we will need two 
lemmas. 



Lemma 4.1 // (^(o;, r) of the form 

ip{u- t) - p ip{y) exp (^-ii^i^^ dy (4.17) 

is a solution of the diffusion equation (4.12), satisfying the initial condition 

(p{uj;t) = ip{uj), 

where (p{i-u) is such that ip{uj, r) is a continuous function of lo for any r > 0, then it 
generates a family of continuous solutions via the transformation 

(f{Lu;T) < — > Aip{2aT-uj;T) e'^°'^°"'~"^\ (4.18) 

where A and a are arbitrary constants. 

Proof. Let us write the function on the right of the relation (4.18) in explicit form: 
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Since 

^ + 2a{aT -io) = - + 2ay, 

It It 

the right-hand side of relation (4.18) is a continuous solution of the diffusion equation 

(4.12) , satisfying the initial condition 

analogously to (4.17). ■ 
The second lemma needed to find the solution of the initial-boundary problem (4.12), 

(4.13) , (4.14) can be stated as follows. 

Lemma 4.2 Consider the function ^p{uj) which verifies to symmetry relation 

ip{u) = ~Lp{2a - uj) e^"^''-'^). (4.19) 

Then, the solution lp{uj]t) of the diffusion equation (4.17), which is continuous with 
respect to lj and with initial condition equal to (p{uj), is vanishing on the line uj = a + aT: 

ip{a + ar; t) = 0, t > 0. 

Proof. Consider the function 

0{UJ- t) = ^{2aT + 2a - W; t) e2a(ar+a-^) ^ (4 20) 

where (/3(w;r) is given by expression (4.17) and ip{iu) obeys to symmetry relation (4.19). 
It follows from (4.16), (4.18) and from the conditions of lemma 4.1, that if{uj; r) satisfies 
the diffusion equation (4.12) and is, for r > 0, a continuous function of the argument 
UJ. Expressions (4.20) and (4.19) ensure that the solution (p{uj;t) satisfies the initial 
condition 

(^(w; T = 0) = ip{2a - uj) e^"'^''-'^^ = ~^{uj). 
This means in turn that 

(^(w;r) = -(^(w;r), 

or in explicit form 

ifiiio; r) = -y.(2ar + 2a - w; r) e2"("-+'^— ). 

In particular 

ip{a + aT;T) = —ip{a + aT;T) => (^(a + ar; r) = 0, r > 0. ■ 

4-5. Solution of the initial-boundary problem 

The solution of the initial-boundary problem (4.12), (4.13), (4.14) is obtained below by 
using the reflection method and the final result is stated in the following theorem. 

Theorem 4.2 The solution of the diffusion equation (4.12), satisfying the initial-boundary 
conditions (4.13), (4.14), is given by 



ip{uj;T)^ J2 e2("-'5)(''-'^)'"'+2("''-/3°)™ X (4.21) 



m— — oo 



g{uj + 2m{b - a); t) - e^''^^^'^ -"'>"'-"'> g{uj + 2m{b - a) - 2a; r) 
© Royal Economic Society 2010 
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where 

9i^;V = ^==cxp - — 
v27rT V 

The proof is given in Appendix A. 2. 
Substituting (4.21) with (4.8) into (4.11), we obtain 



f{h,i,c]K,-f) = g{c- -i) ^ 
where 



^-2(h-lYm^ ~2m(h~'l)(l-K)c 



m— — oo 



1 _ i(h-t}lm-2l(l-(l-K)c) 
1 6 , 

(4.22) 



^.6. Joint pdf of high, low and close values 

Using relations (4.11), (4.15) and (4.22), we obtain the sought joint pdf Q{h,£,c; k,j) 
of the high and low values {H,L} defined by (2.9) of the incomplete bridge Y{t,K,j) 
defined by (2.4), together with the close value C = Ar(l, 7) of the Wiener process X{t, 7) 
with drift given by (2.5). Namely, 

Q{h, £, c; K, 7) = g{c - -f)TZ{h, i; k\c), (4.23) 

where 

cso 

TZ{h,£;K\c)^ m[mV{m{h-i),{l-K)c) + {l-m)V{m{h~i)+i,{l~K)c)] (4.24) 

m — — cxD 

and 

Vih, c) = 4[(c - 2hf - l]e2''('=-''). (4.25) 

Obviously, TZ{h^£;K\c) is the conditional pdf of the high and low values {iJ, L}, under 
the condition that the close value C is equal to a given c. For any k, the conditional 
pdf TZ{h, £] k\c) does not depend on the normalized drift parameter 7. Furthermore, it 
satisfies the normalizing condition 

dh / d£ n{h,£;K\c) = 1. 
^ — 00 

Taking the limit k — > 1 corresponds to the complete bridge, which is an important case 
for our analysis below. Let us thus define the joint pdf limit 

Q{h, £, c; 7) := lim Q{h, £, c; n, 7). 

Expressions (2.5), (4.24) and (4.25) show that it is equal to 

Q{h,£,c;-f) = g{c--f)n{h,£), -00 < c < 00, h > 0, £<0, (4.26) 

where 

00 

n{h,t)= ^ ra[mV{ra{h-£)) + {\-m)V{m{h-l)+l)] (4.27) 

m— — 00 

and 

Vih) = 4(4/1^ - 1) e-'^^\ 

© Royal Economic Society 2010 



Bridge volatility estimators 



13 



Expression (4.26) has a clear probabilistic interpretation. It means that the high and 
low values {H, L} of the complete bridge Y{t, 1, 7) are statistically independent from the 
close value C = X{l,'j) of the original Wiener process with drift. Accordingly, TZ{h,£) 
given by (4.27) reduces to the unconditional joint pdf of the high and low values of the 
complete bridge. 

4-7. Diagrams of the most efficient homogeneous OHLC bridge estimators 

Knowing the joint pdf of the random variables {H, L, C}, one can calculate the auxiliary 
functions gx{0, 0; k, 7) (3.8) needed in the definition of the diagrams (3.10), (3.11) of the 
most efficient homogeneous OHLC volatility bridge estimators. This allows us to derive 
a number of properties of the functions gx{0, (j); k,'^). It follows from (3.8) and (4.23), 
(4.24), that 

gxi9,4>;K,^)= (4.28) 

^ 00 

■]=e^^ 1"^ m [m/A(m(/i — f), c; K, 7) + (1 — to)/a(to(/i — f) + f, c; K, 7)] , 



V27r 
where 

/a(/i, c; k, 7) = y p2+A (^^p _ Y^^) '^{hp^ (1 " n)cp)dp. 

and 

/i = cos 6' cos 0, I — cos s\n (f> , c — sn\9. 
All calculations done, the explicit expression of I\{h, c; k, 7) reads 

3+4 



Ix{h,c;K,'-f) = ( - 



2'2'2ay V2 V2 J \ 2 ' T 2a 

+»r|3 + i)Mr3 + i,|,g)-.„r(2 + i)«(2 + i,|.g 



Here, 



M{a,h,z):=—^ -[ due''''u''~\\-uf-''-\ Re{&} > Rc{a} > 

r(a)r(6 - a) J„ 

is the Kummer fimction (sec Abramowitz M., and A. Stegun. (1964)). We have used the 
following notations 

a = 4h{h - {1 - k)c) + c'^ , 5 = (2/i - (1 - k)c)^ d = 7c. 

In the particular case 7 = 0, we obtain 

Ix{h,c;K):^h{h,c;K,j = 0)= (4.29) 
5+A /3 + X\ (3 + X)[2h - (1 - k)c]2 - {2h - cf - 4:Kch 



[{2h ^ c)^ + AhKc]^ 

Figure 2 shows a 3D plot of the diagram obtained from (3.10), (3.11), (4.28), (4.29) of 
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the most efficient variance bridge estimator, for k = 0.95 and 7 = 0. On the plane (6*, (p) 
depicted bomidary of domain Sk, (3.6). 



1.5 




Figure 2. Diagram of the most efficient variance estimator, for hl — 0.95 and 7 = 



5. COMPARISON OF THE MOST EFFICIENT BRIDGE ESTIMATORS WITH 
THE G&K AND PARK ESTIMATORS 

5.1. Expectation and variance of arbitrary canonical bridge estimators 

In this section, we compare the efficiency of the most efficient homogeneous bridge es- 
timators derived in previous sections with that of the G&K and PARK estimators. We 
thus give the formulas for the expected value and the variance of arbitrary canonical 
homogeneous OHLC bridge estimators defined by (3.4). First, their expected values are 
■A^a(k,7) given by (3.7). 

In general, homogeneous bridge estimators are biased. Thus, one needs a normaliza- 
tion procedure for a practical comparison. We will normalize the homogeneous bridge 
estimators by the one obtained for a zero drift (7 = 0). Thus, for each estimator (3.4), 
we consider its normalized version 

~ _pAVA(e,$) 



Mx{k) 



© Royal Economic Society 2010 



Bridge volatility estimators 15 
where A^a(k) = A^(k,7 = 0). The expected value of the normahzed estimator is 

The variance of the normahzed (at 7 = 0) estimator is 

Var[e.|.,7] = , 

where 

A6.K7)= // i^UO,(l))g2xi0, (I); K,j) cos OdOdcj). 



5.2. Generalized G&K bridge estimator 

We recah that the G&K canonical variance estimator is given by 

ecK - - Lf - k2{C{H + L) - 2HL) - ksC^, (5.1) 
fci = 0.511, fca = 0.019, ks = 0.383. 

The random variables {H,L, C} are the high, low and close values of the Wiener process 
X{t,j) with drift defined by (2.5). In order to compare the efficiencies of the G&K 
estimator and of most efficient bridge estimators, we modify the G&K estimator (5.1) 
by replacing the high, low and close values of the Wiener process X{t,'y) with drift by 
the high, low and close values of the incomplete bridge F(i, ^,7) defined by (2.4). This 
yields 

ecKin) = ki{H- Lf - fc2((l - n)C{H + L) - 2HL) - k^il - nfC^. (5.2) 
The estimator (5.2) can be expressed in a form analogous to (3.4), 

ecK = i?VGK(e,$,K), (5.3) 

with 



V'GK (6*, </>, k) = fci cos^ 9{cos (f> — sin < 



\2 



fc2 



1 — K 

cos^ 9 sin 2(f) — sin 20 (cos (p + sin ( 



— fc3(l — nf sin^ 



To compare the efficiencies of the G&K estimator and of the most efficient bridge 
estimators of arbitrary order A, let us introduce the G&K estimator of order A: 

^~GK,A = TT^^—rT) (0, ^\ (5.4) 

where A^gk,a('*) is given by the following expression 

Mgk,x{k) ^ JJ i:QQ{0,(f),K)gx{0,(j); k) cos 0d0d(j). 

For K ~ and A = 2, the estimator (5.4) reduces to the original G&K estimator defined 
in (5.1). 
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5.3. PARK normalized estimator 
The canonical PARK variance estimator is given by 

We generalize it by the corresponding normalized PARK estimator of order A, 



ep,A = T7-^^p (5-6) 



where 



MO,^) = '"''^41^2'''^'^^ ' Mp,xi^) = 11 i4/\e,cj,)g xie,(b; k) cos 9d9d,^. 



Remark 5.1 Below, we compare the efficiencies of the G&K, of the PARK and of the 
most efficient estimators, and do not discuss the efficiency of the Roger- Satchell estimator. 
Indeed, it follows from our preliminary calculations for k ~ 1 that the Rogers-Satchell 
bridge estimator is significantly less efficient than even the PARK estimator. 

5.4. Comparison of variance estimators 

Figure 3 shows the dependence as a function of the bridge parameter k of the expected 
values of the G&K (5.2) and PARK (5.5) variance estimators, in the case of zero drift 
(7 = 0). One can observe that, for k ^ 0, the G&K and PARK variance estimators are 
biased, so it is convenient to compare their normalized versions (5.4) and (5.6). 




Figure 3. Dependence of the expected values of the G&K (5.2) and PARK (5.5) canonical 
bridge estimators as a function of the parameter k, in the zero drift (7 — 0) case 
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0.2 0.4 



Figure 4. Variances of the most efficient, normalized G&K (5.4) and PARK (5.6) variance 
bridge estimators, as functions of the parameter k, in the case of zero drift (7 = 0) 



Figure 4 plots the numerically calculated dependencies as a function of k, of the vari- 
ances of the most efficient canonical variance bridge estimator, with diagram (3.10), 
(3.11) (A = 2), and of the variances of the G&K and PARK canonical variance estima- 
tors (5.4), (5.6). For k = 0, i.e. in the case of "standard" OHLC estimators, the variances 
of the most efficient and of the G&K estimators are rather close to each other, while the 
variance of the PARK estimator is much larger than the former ones: 

Var[emc,2|K = 0] = 0.2584, Var[eGK,2|K = 0] = 0.2693, Var[gp,2|K 0] = 0.4073. 

In contrast, in the case of an almost complete bridge k e (0.9, 1), the variance of the 
most efficient variance estimator is significantly smaller than the variances of the G&K 
and PARK estimators: 

Var[e„ic,2|K = 1] = 0.1794, Var[gGK,2|K = 1] ^ Var[gp,2|K = I] 0.2. 

Notice that the efficiencies of the G&K and PARK estimators almost coincide for k ~ 1 . 
This is due to the fact that, for k = 1, the G&K variance bridge estimator (5.2) becomes 
close to the PARK estimator: 

gp 2(K = 1) ~ gGK,2(K 1) fci(i?- -Z-)^ + 2A:2i?i, A:i= 0.511, A:2 = 0.0019. 

Simulating 200 realizations of a Wiener process and recording the associated OHLC, 
figure 5 shows the 200 corresponding G&K estimator (5.1), PARK estimator (5.6) (A = 2) 
and most efficient canonical variance bridge estimator, for k = 0.99 and 7 = 0. It 
is visually apparent that the most efficient variance bridge estimator exhibits smaller 
fluctuations and is more efficient than the PARK and G&K bridge variance estimators. 

Remark 5.2 The fact that the OHLC bridge estimators, with k ~ 1, are significantly 
more efficient than the "standard" OHLC estimators, corresponding to k = 0, can be 



© Royal Economic Society 2010 



18 



A. Saichev, D. Sornette, V. Filimonov, F. Corsi 




intuitively explained as follows. It is well-known that the high and low values of a Wiener 
process arc most probably found in the neighborhood of the edges of the observation 
interval. In contrast, by construction of the bridge, its high and low values are in general 
distant from the edges, as illustrated in figure 1. As a result, the high and low of a bridge 
incorporate significantly more information about the behavior of the original stochastic 
process than its own high and low values. 



Remark 5.3 It is noteworthy that the most efficient estimator at 7 = remains more 
efficient than the G&K and PARK estimators as long as 7 remains less than 0.8 (for 
K = 0.95) and similar values for other k's. These condition are not restrictive since 
relevant values of 7 are quite small. Indeed, consider a typical stock with yearly volatility 
a = 0.2 and mean return /i = 0.1. Then, the value of 7 for an estimator calculated at 
the daily scale T ~ 0.004 year is 7 = (^iija^^fT ~ 0.032. For estimators at intra-day 
high-frequencies, for instance for T = 5 minutes = 0.00004 year, we have 7 ~ 0.0032. 
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Figure 6. Variances of the most efficient canonical volatility bridge estimator, G&K (5.4) 
and PARK (5.6) volatility estimators (A = 1), as functions of k (for 7 = 0) 



5.5. Comparison of volatility estimators 

Figure 6 shows the dependencies as a function of k of the variances of the most efficient 
canonical bridge volatility (A = 1) estimator, and the variances of the corresponding 
G&K and PARK volatility estimators (5.4), (5.6). In the case of almost complete bridges 
K e (0.9, 1), the variance of the most efficient volatility estimator is significantly smaller 
than the variances of the analogous G&K and PARK estimators: 

Var[emc,i|« = 1] = 0.0428, Var[gGK,i|«^ = 1] = 0.0473, Var[gp,i|K = 1] = 0.0472. 

Simulating 200 realizations of a Wiener process and recording the associated OHLC, 
figure 7 shows the 200 corresponding G&K volatility (A = 1) estimator (5.4) for k ~ 0, 
the PARK estimator (5.6) (A = 1) and the most efficient canonical bridge volatility 
estimators, for k = 0.99 and 7 = 0. It is visually apparent that the most efficient 
volatility bridge estimator exhibits smaller fluctuations and is more efficient than the 
PARK and G&K bridge variance estimators. 

6. SIMULATED MOST EFFICIENT ESTIMATORS 

The previous sections have derived the most efficient homogeneous bridge estimators, 
whose diagrams (3.10), (3.11), are defined in terms of the function g\{6, cj); k,^) (3.8), 
which depends in turn on the pdf Q{h, £, c; k, 7) (4.23). Notice that relation (3.8) allows 
one to determine the function gx{9, 4>; k, 7) even when the pdf is unknown. The function 
g\{9,(j)]K,"f) can indeed be determined by simulating M 1 times the stochastic pro- 
cess X{t) which describes the statistical properties of the log-price dynamics, and then 
estimate the function g\{6., (j); k,j) by its statistical average. This is particularly conve- 
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Figure 7. Top to bottom: 200 samples of the G&K volatility estimator, for k = 0, PARK 
volatility estimator and most efficient canonical bridge volatility estimator, for k, = 0.99 



nient when theoretical formulas are not available, as occurs when considering stochastic 
processes more complex than the Wiener process with drift. 

To illustrate this possibility of simulating the diagrams associated with the most effi- 
cient estimators, consider the discrete normalized random walk 

f 

X{k) = ^y^e,, k = l,...,K, X(0) = 0, (6.f) 

^=l 

where {e^} is a sequence of iid random variables with zero expectation and unit variance. 
The random walk (6.1) mimics the discrete, tick-by-tick, nature of the log-price stochastic 
process. 

In the limit K — > oo, the random walk X{k) (6.1) converges (even if {ei} are non- 
Gaussian as long as the tail of their pdf is not too heavy) to the Wiener process W{t), so 
that the joint pdf of the random variables {H, L, C} is known theoretically. In contrast, 
in the case of "finite number of ticks" {K < oo), the joint pdf is unknown. Nevertheless, 
one can obtain an approximate expression for the diagram of the most efficient estimator 
by numerical simulation. 

In order to construct the simulated diagram for K = 10; 10^ and 10"^, we divided the 
domain defined in (3.6) in 50 x 50 rectangle bins and, for each K, we generated 
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1.5- 





Figure 8. Synthetic diagram of the most efRcient variance canonical bridge estimator for 
K ~ 10 and k = 1, 7 = 



Table 1. Variances of G&K and simulated most efficient variance estimators. The vari- 
ances of the G&K estimators and of the simulated most efficient variance estimators, 
both for K = and At = 1, are obtained by averaging over N = 10^ simulations of the 
discrete random walk X{k) (6.1). 



K = 10 100 1000 



Var[eGK,2](K = 0) 0.5103 

Var[emc,2](K = 0) 0.4759 

Var[eGK,2](K = 1) 0.4062 

Var[enic,2](K = 1) 0.3373 



0.3272 0.2858 0.2693 

0.3130 0.2755 0.2584 

0.2434 0.2125 0.1996 

0.2151 0.1896 0.1794 



M — 10^ simulations of the random walk X{k) (6.1) with Gaussian summands {ci}- 
We then calculated the function g\{0, 4>; k) (for 7 = 0) using the approximate statistical 
relation 

1 

gx{e,<P; k) cos9d9d(t)~ — ^ R^Js{<dm,'^m). 

771 — 1 
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Figure 9. Top to bottom: numerical samples of the G&K estimator, for k = and 
K = 1, and the numerical samples of the simulated most efficient canonical bridge variance 
estimator, for k = 1, obtained for discrete random walk X{k) (6.1) with K = 10 



The set {8™, 'J*m, Rm} are samples of the random variables (3.3) obtained for the m-th 
simulation, is the indicator of the set 

5= {6,6 + de) X {(j),(j) + d(j)), 

and the summation is performed over M simulations of the random walk (6.1). The 
histograms of the function g\{9, </>; k) thus obtained is then substituted into the diagram 
function (3.10), (3.11) to produce its 2D linear interpolation. 

Figure 8 presents the 3D plot of the synthetic diagram of the most efficient variance 
estimator, obtained by statistical averaging for K = 10 and k = 1. 7 = 0. Notwithstand- 
ing the visible fluctuations, table 1 shows that this level of numerical approximation is 
sufficient to obtain significantly better efficient estimators than for instance, G&K esti- 
mator. Table 1 gives the variances of the canonical variance bridge estimators. The values 
shown in table 1 have been obtained by statistical averaging over iV = 10^ simulations 
of the random walk (6.1). 

Simulating 200 realizations of a Wiener process (6.1) for K = 10, figure 9 shows the 
simulated most efficient canonical variance bridge estimator, for /t = 1, and the samples 
of the G&K estimators, for k = and k = 1. It is clear that the simulated most efficient 
bridge variance estimator is significantly more efficient than the G&K estimator. 
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7. CONCLUSIONS 

In this paper, we have pursued the development of a comprehensive theory of homoge- 
neous volatihty estimators of arbitrary stochastic processes. Our focus has been to derive 
OHLC (opcn-high-low-closc) log-prices bridge volatility estimators, which can span time 
intervals extending from seconds to years. The main tool of our theory is the parsimonious 
encoding of all the information contained in the OHLC in the form of general "diagrams" 
associated with the joint distributions of the high-minus-open, low- minus-open and close- 
minus-open values of the original log-price process and its bridge. The diagrams can be 
tailored to yield the most efficient estimators associated to any statistical properties of 
the underlying log-price stochastic process. 

Previous works have developed variance estimators which are quadratic functions of the 
OHLC. Our main contribution is to stress the remarkable fact that quadratic estimators 
are only particular cases of general homogeneous estimators. Our theory constructs the 
tools to find most efficient homogenous estimators which, by construction, are always 
more efficient than the most efficient quadratic estimators. Perhaps paradoxically, it 
turns out that the search for the most efficient quadratic estimators is more tedious than 
that of the more efficient homogeneous estimators. Another advantage of homogeneous 
estimators is that they give the possibility to develop efficient volatility in addition to 
variance estimators, while quadratic estimators are specialized to variance estimators. 

Our theory opens several interesting developments. First, the determination of the key 
functions g\{9^ (j);"f), defining the above diagrams, provides the tools to develop efficient 
bridge volatility estimators for arbitrary non-Gaussian log-price processes, including the 
presence of micro-structure as in tick-by-tick price series. Our methods should lead to the 
development of effective algorithms for low- and high-frequency OHLC volatility bridge 
estimators, that can be applied in practice to any kind of financial markets. 
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A. APPENDIX 
A.l. Proof of Theorem 3.1 

Eor given values of the parameters n and 7 = 70, the variance of the unbiased canonical 
estimator, with diagram (3.9), is equal to 



Var[eA;K,7o] = 




(A.l) 



Using the Schwarz inequality 
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for arbitrary locally integrable real- valued functions A{9, (p) and B{d, (j)), we take 



Ai9, 0) = Gie, 0)V52A(e,0;«:,7o)cos0, 



cos ( 



g2A(6',0;K,7o)' 



and obtain 




G{9,(j))gx{e,(j)] K,j) cos 




G\e,cj,)g2x{eA;^.lo) cos OdOdcj, II M^MlJhJ^ COS 

5'2a(6',(/>;k,7o) 

It follows from (A.l) and from the above inequality that the variance of any canonical 
homogeneous volatility estimator of order A satisfies the inequality 

Var[eA;K,7o] ^Fa(k,7o), ^4(^^,7)^ / - 1, (A.2) 

where £x{k, 7) is defined by expression (3.11). It follows from (A.l), (A.2) and (3.11) that 
the variance of the canonical volatility estimator of order A reaches its minimal value for 
a given 7 = 70 and k. if G{0, cf)) is given by the left equality of (3.11). ■ 

A.2. Proof of Theorem 4-- 2 

It is convenient to replace the initial condition (4.13) by the more general one 

(p{uj; T = 0) = Lp{uj), ijj e (a, h). (A. 3) 

At the end of proof, we obtain formula (4.21) by taking '^{ijj) — 5{uj). 

The idea of the proof consists in redefining the function (p{uj) in (A. 3) outside the 
interval lj S (a, 6) in such a way that the solution of equation (4.12), supplemented by 
the initial condition 

(.p{uj]T ^ Q) = (f{u), ljg(— 00,00), (A. 4) 

satisfies the absorbing boundary conditions (4.14). In other words, it should be equal to 
zero on the lines lj ^ a + ar, lj = h -\- Pt, r > 0. Let us define the auxiliary function 

(Po{lj) ^ ip{Lj)l(^a,b){^^), t^e (-00,00), (A.5) 

where Ie{x) is the indicator of the set E. 

It follows from lemma 4.2 that the solution of the diffusion equation, supplemented 
by the initial condition (A. 4), satisfies the boundary conditions (4.14) if (p{lj) obeys to 
symmetry relations 

ip{Lj) = -ip{2a - Lo) e^"^''-'^\ ifiio) = -ip{2b - Lo) e^^'^''-'^\ (A.6) 

Using the first of these two equalities and definition (A.5) of the function ipo{Lo), let us 
redefine ip{lj) onto the interval oj G (2a — b, b) as follows: 

ip{Lj) ^ ifi'^iuj), we (2a -6, 6), /(w) = (po('^)-<Po(2a-w) e2"('^-'^). 
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Then, the equahties (A. 6) provide the "quasiperiodic" relation 

^{u;) ^ ^{u; + 2{b - a)) e'(P-^){^+b-a)+2(<.b-pa) ^ 

which yields 

oo 

m— — oo 

where 

= /(CJ + 2(6 - a)m) e2(/3-'^)(-'+™A)m+2(a6-/3a)m_ 

Substituting </j(w) given by (A. 7) with (A. 8) into (4.17), we obtain the sought solution 
of the initial-boundary problem (4.12), (4.13), (4.14). In particular, using ^Pvi{lo) = <5(a;), 
that is 

^"(w) ^°(cj) = 5(w)-e-2"° (5(w-2a), 

we obtain the solution (4.21). ■ 
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